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ABSTRACT 

We conduct a comprehensive investigation of the detectability of the first stars and their enrichment 
signatures in galaxy clusters. As the initial mass function (IMF) of these Population III (PopIII) stars 
is unknown, and likely to be biased to high masses, we base our study on analytical models that param- 
eterize these uncertainties and allow us to make general statements. We show that the mean metallicity 
of outflows from PopIII objects containing these stars is well above the critical transition metallicity 
(Z CT ~ 10 -4 Zq) that marks the formation of normal stars. Thus the fraction of PopIII objects formed 
as a function of redshift is heavily dependent on the distribution of metals and fairly independent of the 
mean metallicity of the universe, or the precise value of Z CI . 

Using an analytic model of inhomogeneous structure formation, we study the evolution of PopIII 
objects as a function of the star formation efficiency, IMF, and efficiency of outflow generation. For 
all models, PopIII objects tend to be in the 10 6 5 M Q - 1O 7 -°M mass range, just large enough to 
cool within a Hubble time, but small enough that they are not clustered near areas of previous star 
formation. Although the mean metallicity exceeds Z CI at z ~ 15 in all models, the peak of PopIII star 
formation occurs at z ~ 10 and such stars continue to form well into the observable range. We discuss 
the observational properties of these objects, some of which may have already been detected in ongoing 
surveys of high-redshift Lyman-a emitters. 

Finally, we combine our PopIII distributions with the yield models of Heger and Woosley (2002) to 
study their impact on the intracluster medium (ICM) in galaxy clusters. We find that PopIII stars can 
contribute no more than 20% of the iron observed in the ICM, but if they form with characteristic masses 
~ 200 — 260M Q , their peculiar elemental yields help to reconcile theoretical models with the observed 
Fe and Si/Fe abundances. However, these stars tend to overproduce S/Fe and can account only for the 
O/Fe ratio in the inner regions of poor clusters. Additionally, the associated SN heating falls far short 
of the observed level of ~ IkeV per ICM gas particle. Thus the properties of the first objects may be 
best determined by direct observation. 



1. INTRODUCTION 

The properties of the first cosmic stars born out of a pri- 
mordial, metal-free gas may have been quite different from 
those of later generations formed from polluted material. 
Studies of the gravitational collapse of the 1O 6_7 M0 pre- 
galactic units that harbor the first stars suggest that these 
Population III (PopIII) objects fragmented into clumps 
with typical masses of 1O 3_4 M0 (Abel et al. 1998; Bromm, 
Coppi, & Larson 1999, 2002; Abel, Bryan & Norman 
2000; Bromm et al. 2001) that did not subdivide further 
(Omukai & Nishi 1998; Nakamura & Umemura 1999; Ri- 
pamonti et al. 2002). As the strong accretion rates (up to 
0.1 M©yr _1 ) onto the central protostellar cores in these 
clumps could not have been arrested by radiation pressure, 
bipolar outflows, rotation or by any other known mecha- 
nism (Ferrara 2001), they are likely to have evolved into 
very massive stars. Such massive star formation mode 
would have continued until the gas metallicity passed a 
critical threshold of Z cr 10~ 4 Zq (Omukai 2000; Bromm 
et al. 2001; Schneider et al. 2002), which enabled both the 
formation of lower mass clumps and some of the above 
mechanisms to stop accretion. 

This situation leads to the so-called star formation co- 



nundrum pointed out in Schneider et al. (2002). If PopIII 
stars were indeed very massive M* > IOOMq, stellar evo- 
lution models predict that most of these objects would 
have evolved into black holes, which lock in their nu- 
cleosynthetic products. This in turn would prevent the 
metal enrichment of the surrounding gas and force the 
top-biased star formation mode to continue indefinitely. 
It is only in a relatively narrow mass window in fact, 
A 77 = 140M Q < M* < 260M Q , that ejection of heavy el- 
ements can occur, powered by so called pair-instability su- 
pernovae in which a pair-production instability ignites ex- 
plosive nucleosynthesis. Thus these supernovae may be es- 
sential to increase the intergalactic medium (ICM) metal- 
licity above Z cr and initiate the transition from top-biased 
star formation to the formation of "normal" (PopII/PopI) 
stars with a more typical initial mass function (IMF). 

While this transition must have occurred at high redshift 
in the bulk of cosmic star formation sites (as suggested by 
the predicted diffuse neutrino flux associated to the black 
hole collapse of the first stars, Schneider, Guetta, & Fer- 
rara 2002, and by the interpretation of the near infrared 
background, Salvaterra & Ferrara 2003) massive PopIII 
stars are likely to have left some imprints that are observ- 
able today. One point that might have important observa- 
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tional consequences is the fact that cosmic metal enrich- 
ment has proceeded very inhomogeneously (Scannapieco, 
Ferrara, & Madau 2002; Thacker, Scannapieco, & Davis 
2002; Furlanetto & Loeb 2002), with regions close to star 
formation sites rapidly becoming metal-polluted and over- 
shooting Z cr , and others remaining essentially metal-free. 
Thus, the two modes of star formation, PopIII and nor- 
mal, must have been active at the same time and possibly 
down to relatively low redshifts, opening up the possibility 
of detecting PopIII stars. 

This is particularly important as metal-free stars are 
yet to be detected. This could be attributed to a sta- 
tistical paucity of metal-free stars or may be due to 
their absence entirely, hinting at short-lived, massive first 
stars. At the same time, theoretical studies are making 
rapid progress in understanding the radiative properties 
(Tumlinson, Giroux, & Shull 2001; Bromm, Kudritzki, 
& Loeb 2001) and spectral energy distributions of such 
stars (Schaerer & Pello 2002; Oh, Haiman, & Rees 2001), 
and their interactions with the ambient medium (Ciardi 
& Ferrara 2001). Thus, it might well be that the best 
approach to detect the first stellar clusters is to develop 
careful observation strategies of high-redshift PopIII host 
candidates. This is the first aspect that we investigate in 
detail in this paper. 

However, direct observations are not the only way to 
learn about the first cosmic stars. Abundance patterns 
measured in different cosmic objects may already show 
clear imprints of such early star formation. A prime ex- 
ample are observations of metal-poor halo stars. Although 
these stars cannot be the massive PopIII stars themselves, 
they might be old enough to retain some record of the 
metals ejected by the first stellar generation (Cayrel 1996; 
Weiss et al. 2000). Indeed, high-resolution spectroscopy of 
metal poor-halo stars has uncovered a number of unusual 
abundance patterns in the iron peak elements in stars with 

[Fc/H] 2.5 (McWilliam et al. 1995; Ryan ct al. 1996). 

In these stars the mean values of [Cr/Fe] and [Mn/Fe] de- 
crease toward smaller metallicities, [Co/Fe] increases with 
metallicity, and r- and s-process elements are nearly ab- 
sent: trends that cannot be interpreted with Type II or 
Type la SNa yields (Oh et al. 2001; Umeda & Nomoto 
2002). Very recently, the Hamburg/ESO objective prism 
survey (HES) team has reported the discovery of a star 
HE0107-5240 with a mass of 0.8M Q and an iron abun- 
dance of [Fe/H] = -5.3 ± 0.2 (Christlieb et al. 2002). 
While this iron value is uniquely low, the overall metal- 
licity of this star is within the estimated range for Z cr 
(1O- 6 Z < Z cr < 10- 4 Z Q , Schneider et al. 2002, 2003) 
suggesting that the identification of more stars with iron 
abundances less than [Fe/H]< —5 will be extremely im- 
portant to elucidate the conditions that enable low-mass 
star formation. 

A second indirect approach is to search for PopIII en- 
richment in the intcrgalactic medium itself. Tremendous 
progress has been made recently in the determination of 
the metal content of the IGM both observationally (El- 
lison et al. 2000; Schaye et al. 2000; Songaila 2001) and 
theoretically (Gnedin & Ostriker 1997; Ferrara, Pettini, 
& Shchckinov 2000; Ccn & Bryan 2001; Scannapieco & 
Broadhurst 2001; Madau, Ferrara, & Rees 2001; Scanna- 
pieco, Ferrara, & Madau 2002; Thacker, Scannapieco, & 
Davis 2002; Furlanetto & Loeb 2002). However, abun- 



dance patterns in the low-column density Lya forest can 
be obtained only for a few elements and are complicated 
by uncertain ionization corrections. 

In this work, we examine a third possibility for searching 
for signatures of PopIII star formation, by studying the in- 
tracluster medium (ICM). Metal abundance experiments 
are easier in the ICM than in the more diffuse IGM due not 
only to higher gas densities, but to the lower redshift of the 
targets, the simpler physics of the collisionally ionized gas, 
and the larger statistics. Such an approach parallels simi- 
lar efforts by Renzini (1999) and Loewenstein (2001), but 
it is specifically designed to test for PopIII star formation. 

The properties of ICM may also provide information 
as to the energetics of SN 77 , whose high energy output 
is suggestive of well-know entropy-floor problem (Loewen- 
stein 2000; Lloyd-Davies, Ponman, & Cannon 2000; Tozzi 
& Norman 2001; Borgani et al. 2001). Here a wide range 
of observations have shown that ICM is much more diffuse 
and extended than predicted by simple virialization mod- 
els, resulting in an X-ray luminosity-temperature relation 
that is widely discrepant with simple predictions (Kaiser 
1991). This discrepancy appears to be due to a baseline 
entropy of 55— 110 /i -1 / 3 kcV cm 2 , where this uncertainty 
is due to the unknown extent to which shocks have been 
suppressed in low-mass systems. This excess appears to 
be distributed uniformly with radius outside the central 
cooling regions and provides an additional constraint on 
our models of PopIII formation. 

Our aim is hence twofold. We first investigate the spa- 
tial distribution and redshift evolution of galaxies hosting 
PopIII stars. This study enables us to predict a number 
of observationally related quantities (detection probabili- 
ties, number counts, equivalent width distributions) that 
we combine into a search strategy to find these systems. 
Secondly, we turn to the effects of these metal- free stars on 
the metallicity patterns and energetics of galaxy clusters. 

Throughout this paper, we restrict our attention to a 
cosmological model in which h = 0.65, CIm — 0.3, Qa = 
0.7, n b = 0.05, <r$ = 0.87, T = 0.18, and n = 1, where 
h is the Hubble constant in units of 100 km s _1 Mpc -1 , 
Qm, and fib are the total matter, vacuum, and bary- 
onic densities in units of the critical density, as is the mass 
variance of linear fluctuations on the 8ft. _1 Mpc scale, T is 
the CDM shape parameter, and n is the tilt of the pri- 
mordial power spectrum. The choice of these parameters 
is based mainly on measurements of Cosmic Microwave 
Background anisotropies (eg. Balbi et al. 2000; Netterficld 
et al. 2002; Pryke et al. 2002) and the abundance of galaxy 
clusters (eg. Viana & Liddle 1996). 

The structure of this work is as follows. In §2 we derive 
metal yields and kinetic energy input from pair-creation 
and Type II SNe, and in §3 we use these quantities to 
construct a model of metal dispersal by SN-driven out- 
flows. In §4 we present an analytic model of inhomo- 
geneous structure formation which we use in §5 to esti- 
mate the distribution of cosmic objects hosting PopIII or 
PopII/I stars, accounting for feedback effects due to out- 
flows pre-enriching neighboring halos. In §6 we apply these 
results to construct the Lyman-a flux and redshift distri- 
bution of PopIII objects and discuss strategies for their 
direct detection. In §7 we turn our attention to the intra- 
clustcr medium of galaxy clusters, and relate the thermal 
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and enrichment properties of this gas to the properties of 
PopIII objects. Conclusions are given in §8. 

2. METAL YIELDS AND EXPLOSION ENERGIES 

In order to quantify the properties of PopIII and PopII/I 
objects, we consider the metal yields and explosion en- 
ergies of pair-creation (SN 77 ) and Type II SNe (SNII), 
respectively. In the PopII/I case we adopt the results 
of Woosley & Weaver (1995) who compute these quan- 
tities for SNII as a function of progenitor mass in the 
range 12M Q < M* < 40M Q and initial metallicities 
Z = (0, 1CT 4 , l(T 2 , 1CT 1 , 1) Z Q . The corresponding quan- 
tities for SN 77 are taken from Heger & Woosley (2002) 
assuming progenitors in the mass range A 77 = 140M Q < 
Mi, < 26OM . 

2.1. PopII/I Objects 

The total heavy-element mass and energy produced 
by a given object depend on the assumed IMF. We de- 
fine PopII/I objects as those that host predominantly 
PopII/PopI stars, i.e. objects with initial metallicity Z > 
Z CT = 10~ 4 Z Q . For these, we assume a canonical Salpeter 
IMF, $(M) oc M-( 1+x \ with x = 1.35, and lower (upper) 
mass limit M t = 0.1M Q (M u = 100M Q ), normalized so 
that, 

dMM$(M) = 1. (1) 
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The IMF-averaged SNII metal yield of a given heavy ele- 
ment i (in solar masses) is then 



CL dM<&(M)Mi 



(2) 



where M, is the total mass of element i ejected by a pro- 
genitor with mass M. The mass range of SNII progenitors 
is usually assumed to be (8 — 1OO)M . However, above 
Mbh = 50±10M Q stars form black holes without ejecting 
heavy elements into the surrounding medium (Tsujimoto 
et al. 1995), and from (8 — 11)M© the pre-supernova evo- 
lution of stars is uncertain, resulting in tabulated yields 
only between 12M and 40M© (Woosley & Weaver 1995). 
To be consistent with previous investigations (Gibson, 
Loewenstein, & Mushotzky 1997; Tsujimoto et al. 1995), 
in our reference model (SNII-C in Table 1) we consider 
SNII progenitor masses M S nii = 10M Q < M* < M BH = 
50Mq range, linearly extrapolating from the lowest and 
highest mass grid values. Finally, SNII yields depend on 
the initial metallicity of the progenitor star. Here we as- 
sume for simplicity that SNII progenitors form out of a gas 
with Z > 10~ 2 Z Q . This choice is motivated by the typical 
levels of pre-enrichment from PopIII stars (see Section 3). 
Moreover, the predicted SNII yields with initial metallic- 
ity in the range 10~ 4 Z Q < Z < 10~ 2 Z & are largely inde- 
pendent of the initial metallicity of the star (Woosley & 
Weaver 1995). Different combinations of progenitor mod- 
els relevant for the present analysis are illustrated in Table 
1. The last entry in each row is the number of SNII pro- 
genitors per unit stellar mass formed, 
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Finally, we assume that each SNII releases £kin = 1.2 x 
10 51 erg, independent of the progenitor mass (Woosley & 
Weaver 1995). 

Table 3 shows the IMF-averaged yields for some ele- 
ments relevant to our investigation as well as the total mass 
of metals released in different SNII models. In the same 
Table, we also show the corresponding yields for Type 
la SNe (SNIa), whose values are mass- independent and 
have been taken from Gibson, Loewenstein, & Mushotzky 
(1997). 

2.2. PopIII Objects 

The detailed shape of the PopIII IMF is highly uncer- 
tain, as these stars have not been directly observed. In 
this investigation we adopt a model that is based on nu- 
merical and semi-analytical studies of primordial gas cool- 
ing and fragmentation (Omukai 2000; Bromm et al. 2001; 
Nakamura & Umemura 2001; Schneider et al. 2002; Abel, 
Bryan, & Norman 2002), which show that these processes 
mainly depend on molecular hydrogen physics. In partic- 
ular, the end products of fragmentation are determined by 
the temperature and density at which molecular hydrogen 
levels start to be populated according to LTE, significantly 
decreasing the cooling rate. The minimum fragment mass 
is then comparable to the Jeans mass at these conditions, 
which is ~ 1O 3 M0, although minor fragmentation may 
also occur during gravitational contraction, due to the oc- 
casional enhancement of the cooling rate by molecular hy- 
drogen three-body formation. 

Ultimately, in these conditions the stellar mass is deter- 
mined by the efficiency of clump mass accretion onto the 
central protostcllar core. For a gas of primordial composi- 
tion, accretion is expected to be very efficient due to the 
low opacity and high temperature of the gas. Given the 
uncertainties is these models, it is interesting to explore 
different values for the characteristic stellar mass Mc and 
for ere, the dispersion around this characteristic mass. In 
the following we assume that PopIII stars are formed ac- 
cording to a mass distribution given by, 



$(M)MdM 



1 



e -(M-M c ? /2*l dM ^ (4) 



where 100M o < M c < 1000M Q and <j™' m < <j c < a% ax . 
A similar shape for the initial PopIII IMF has been also 
suggested by Nakamura & Umemura (2001), with a second 
Gaussian peak centered around a much smaller character- 
istic mass of IMq, leading to a bimodal IMF. 

The minimum dispersion is taken to be <7™ ln = O.IMq 
and the maximum value is set so that less than 1% of 
stars form below a threshold mass of Mbh ~ 50M Q , i.e. 
a™ 3 -* = (M c - 5OM )/3, such that they will not leave be- 
hind small stars or metal-free compact objects observable 
today. With this choice, only progenitors with mass in the 
range A 77 are effective in metal enrichment. The IMF- 
averaged yields, Y^ 77 , and the number of progenitors per 
unit stellar mass formed, A/" 77 can be obtained from ex- 
pressions equivalent to eqs. (2) and (3) with the assumed 
IMF and metal yields from Heger & Woosley (2002). For 
each pair of Mq and ac values, we also compute the pa- 
rameter 

f A dMM$(M) 
777 " f °° dMM^(M) ' 1 j 
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Table 1 

The different progenitor models for SNII. The masses arc in Mq and M 11 is in units of Mq -1 . 



SN Model 


Msnii 


M BH 


Z/Z Q 


N n 


SNII-A 


12 


40 


io- 2 


0.00343 


SNII-B 


12 


40 


1 


0.00343 


SNII-C 


10 


50 


io- 2 


0.00484 


SNII-D 


10 


50 


1 


0.00484 



Table 2 

The first three entries are in Mq, A/" 77 is in units Mq" 1 and ftin is m units of 10 51 erg. The last entry 
represents the kinetic energy per unit gas mass, = fl 11 f m N" n '££?J in units of 10 51 erg Mq" 1 computed assuming fl 11 = 0.1 and 

f w = 0.1 (see text). 



SN Model 


M c 


G c 


a C 


A/" 77 


fn 


c kin 


clll 


SN 77 -A 


100 


10 


17 


< 6 x IO" 5 


[0.03- 8] IO" 3 


6 - 7.5 


< 4 x IO" 6 


SN 77 -B 


200 


20 


50 


[5-4]10~ 3 


0.99-0.8 


38.4- 37.5 


[2- 1.5]10" 3 


SN 77 -C 


260 


26 


70 


[2.1 - 2.2]10~ 3 


0.5-0.4 


67.1 -46.4 


[1.4- 1]10" 3 


SN 77 -D 


300 


30 


83 


[0.4- 1.4]10- 3 


0.09-0.3 


70.5 - 48.2 


[2.6- 6.6]10~ 4 


SN 77 -E 


500 


50 


149 


< 2.2 x 10" 4 


< 0.05 


70.6 - 47.8 


< io- 4 


SN 77 -F 


1000 


100 


314 


< 3 x IO" 5 


< 6.3 x 10~ 3 


69.9 - 43.6 


1.3 x IO -5 



Table 3 

The IMF-averaged metal yields (in Mq) for SNII, SNIa, and SN. 



SN Model 


Yo 




Y S 


Y Fc 


^rnet 


SNII-A 


1.43 


0.133 


0.064 


0.136 


2.03 


SNII-B 


1.56 


0.165 


0.078 


0.115 


2.23 


SNII-C 


1.11 


0.1 


0.047 


0.126 


1.6 


SNII-D 


1.19 


0.124 


0.061 


0.11 


1.76 


SNIa 


0.148 


0.158 


0.086 


0.744 


1.23 


SN 77 -A 


48.4-47.5 


2.07-4.37 


0.559-13.2 


[0.29 - 13.3]10" 2 


59.3 - 61.3 


SN 77 -B 


44.2-43.5 


20.9-18.9 


8.77-7.89 


5.81-7.41 


89 - 86.8 


SN 77 -C 


37.4-41.8 


24.3-20.8 


11.2-9.01 


22.1-11.6 


104 - 92.4 


SN 77 -D 


35-41.4 


23.6-21.1 


11.0-9.19 


24.9-12.5 


104 - 93.4 


SN 77 -E 


32.8-41.3 


22.5-20.8 


10.5-9.07 


25.7-12.4 


100 - 92.9 


SN 77 -F 


34.3-42.1 


23.1-19.8 


10.8-8.51 


24.8-10.6 


102 - 90.2 
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i. e. the fraction of PopIII stars ending in SN 77 . Similarly, 
the IMF-averaged specific kinetic energy is 



r;n _ J Ayi dM$(M)E kin 



cm 

c kin = 



LdM^(M) 



(6) 



The various models are described in Table 2 and the cor- 
responding yields are given in Table 3. 

3. GALAXY OUTFLOWS AND METAL DISPERSAL 

Having constrained the yields and energetics of individ- 
ual starbursts, we now relate them to the properties of the 
resulting SN driven outflows. Here we simply model out- 
flows from both PopIII and PopII/I objects as spherical 
shells expanding into the Hubble flow, as in the formal- 
ism described by Ostriker & McKee (1988) and Tegmark, 
Silk, & Evrard (1993). These shells are driven only by 
the internal hot gas pressure and decelerated by inertia 
due to accreting material and gravitational drag while es- 
caping from the host. As the properties of PopIII objects 
are quite uncertain, the inclusion of more detailed physical 
processes such as cooling, stochastic variations in the star 
formation rate (SFR), and external pressure are not war- 
ranted for our purposes here (eg., Scannapieco & Broad- 
hurst 2001; Madau, Ferrara, & Rees 2001; Scannapieco, 
Fcrrara, & Madau 2002). In this case the relevant evolu- 
tionary equations become 



Rs 
E b 



pR s R s 



HR a f - Q 



H R s 



M " 



L(t) - ^R 2 s R s P b , 



(7) 



where the overdots represent time derivatives, the sub- 
scripts s and b indicate shell and bubble quantities respec- 
tively, R s is physical radius of the shell, E b is the internal 
energy of the hot bubble gas, Pi, is the pressure of this gas, 
and p is the mean IGM background density. Finally we 
assume adiabatic expansion with index 7 = 5/3 such that 
P b = E b /2nR 3 s . 

The evolution of the bubble is completely determined by 
the mechanical luminosity evolution assigned, L(t). Again 
for simplicity, we take both PopIII and PopII/I objects to 
undergo a starburst phase of tsN = 5 x 10 7 years, but with 
different energy- input prescriptions. We approximately 
account for the gravitational potential of the host galaxy 
by subtracting the value of GM 2 (0(,/fiM)/^vir, from the 
total wind energy, where M is the total mass (dark + 
baryonic) of the object and r vlT is its virial radius. Taking 
the standard collapse overdensity value of 178 this gives a 
mechanical luminosity of 

L(t) = 16OL M 6 6(i S N - t) [fl 1 ' 111 'AAinJV"' 77 

-5 xlO- 12 (M b M7 M M) 2/3 (l + *<•)] , (8) 

where /+ is fraction of gas converted into stars, f w is the 
fraction of the SN kinetic energy that is channeled into 
the galaxy outflow, M b is the baryonic mass of the galaxy 
in units of solar mass, and z c is the collapse redshift of 
the object. Taking a typical estimate for PopII objects of 
f" fwZkinN 11 = 10- 4 '(in units of 10 51 ergsM©" 1 ), we see 
that the gravitational contribution becomes important for 
these objects when M b hn M /^ b £ 10 n (l + z c )~ 3/2 M . 



The only difference in the wind evolution of PopIII 
and PopII/I arises from the product, fl 1 ' 111 fwE-kinN 11,111 , 
which we define as the "energy input per unit gas mass" 
£g I,IH , which has units of lO^ergM" 1 . In the PopII/I 
case, we take // 7 to be 0.1, which gives good agreement 
with the observed high redshift star formation rates and 
abundances of metals measured in high-redshift quasar 
absorption line systems (Thacker, Scannapieco, & Davis 
2002; Scannapieco, Ferrara, & Madau 2002). Also in this 
case we constrain f w by combining the overall efficiency of 
30% derived for the 2 x 10 8 M Q object simulated by Mori, 
Fcrrara, & Madau (2002) with the mass scaling derived 
in Ferrara, Pettini, & Shchckinov (2000), which was ob- 
tained by determining the fraction of starburst sites that 
can produce a blow-out in a galaxy of a given mass. This 
gives f w (M) = 0M B (M)/S B (M = 2 x 1O 8 M ) where 



6 B (M) 



1.0 

1.0- 0.165 HNf 1 ) 
[1.0-0.165 ln(100)] 100 A> t " 



N t < 1 

1 < Nt < 100 , 
100 < N t 

(9) 

and N t = 1.7 x 10 7 (tt b /ttM)M/M Q is a dimensionless 
parameter that scales according to the overall number of 
SNe produced in a starburst, divided by the star formation 
efficiency, fl 1 . 

In the PopIII case on the other hand, there arc no di- 
rect constraints on either f^ 11 or the wind efficiency. For 
these objects we allow these parameters to be free, varying 
Eg 11 over a large range as discussed below. Finally, when 
outflows slow down to the point that they are no longer 
supersonic, our approximations break down, and the shell 
is possibly fragmented by random motions. At this point 
we let the bubble expand with the Hubble flow. 

Using this model, we can compute the comoving ra- 
dius at a redshift z r of an outflow from a PopIII 
source of total mass M and formation redshift z Sl 
rin(M,z s ,z r ) = rlflh/QMMSg 11 , z s , z r ], and the equiv- 
alent quantity for PopII/I objects, rn(M,z s ,z r ) = 
r^h/QuMEg 1 [M), z s , z r \. In the upper panels of Fig. 1, 
we plot nn and rjj as functions of z r for sources of various 
mass scales, choosing two representative E^ 11 values drawn 
from Table 2. In the left panels we consider lower redshift 
(z s = 10) sources with masses 10 7 M Q to 1O 9 M in a model 
with Eg 11 = 10~ 4 ' , while in the right panels we consider 
more-energetic (E^ 11 = 10~ 2 - 5 ), earlier (z s — 20) sources 
with masses from 1O 6 M to 1O 8 M . In the PopII/I case, 
we take f+EunN 11 — 4.84 x 10~ 4 , consistent with models 
SNII-C and SNII-D, with a 10% star formation efficiency 
and a supernova energy of w 10 51 ergs. Here we see that 
each source ejects material to beyond the radius of the 
perturbation from which it was formed, and therefore this 
material is easily able to affect the formation of neighbor- 
ing objects. 

In the lower panels of Fig. 1, we compute the mean 
metallicity within the expanding bubbles, for the same 
range of parameters as in the upper curves. For simplicity, 
we assume that in the PopIII case, the mass released in 
metals is equal to twice the kinetic energy content, that is 
Ymet/Mg) w 2£kin/10 51 erg, which is a reasonable approxi- 
mation for the models in Tables 2 and 3. Throughout this 
paper, we adopt a simple model that assumes that half of 
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the metals are ejected into the outflow in both the PopII/I 
and the PopIII cases. Furthermore, we adopt a single 
metallicity for this material, ignoring what is likely to be 
a strong gradient from the bubble to the shell. Finally, we 
do not attempt to reproduce the detailed interaction be- 
tween SNe and the inhomogcneous ISM of the galaxy, an 
issue studied in detail in Nakasato and Shigeyama (2000). 

For the purposes of Fig. 1, we take f w = 1.0 in all PopIII 
objects, which translates into a minimum metallicity in 
each outflow, as this gives the largest values of S^ 11 for a 
fixed £kin value. In the PopII case, metal masses are com- 
puted as per model SNII-C. The final metallicity of the 
bubble scales roughly as M 2 / 5 as the total mass in metals 
in the bubble is oc M while its final comoving volume is 
oc M 3 / 5 , due to the r oc E x / b Sedov- Taylor scaling. Thus 
increasing the object mass by a factor of a hundred in- 
creases the final wind metallicity by almost an order of 
magnitude, an effect that is enhanced further in PopII 
objects by the assumed inefficiency of large outflows, as 
parameterized by eq. (9). However, even in the smallest 
PopIII objects, with the smallest assumed values of E^ 11 
and metal masses, the mean metallicity of the resulting 
bubble is well above the transition metallicity Z CT . As the 
parameters in this figure were chosen to provide a strong 
lower limit on the metallicity of these systems, a clear im- 
plication is that all enriched regions of the IGM lie above 
the critical value necessary for the formation of PopII/I 
stars. 

Hence, the switch between primordial and recent star 
formation was determined primarily by the spatial distri- 
bution of forming sources and metals in the universe. This 
is an essential point in understanding the properties of pri- 
mordial enrichment, which is independent of the details of 
our models, as the lowest levels of Z w ind derived from our 
estimates are several times higher than the most likely 
values for Z CT . The fraction of PopIII objects formed as a 
function of redshift is intimately tied up with the evolu- 
tion of the outflows themselves and fairly independent of 
the mean metallicity of the universe, or the precise value 
of the critical transition metallicity. Thus, the evolution 
of PopIII objects can only be properly tracked from an ap- 
proach that captures the spatial distribution of collapsing 
halos, as such objects are able to form in pristine areas over 
a large range of redshifts. Conversely, primordial objects 
are able to enrich selected areas of the IGM to levels well 
above Z cr and may have made a significant contribution 
to the metals found in the biased areas in which galaxies 
formed most vigorously. 

4. SPATIAL DISTRIBUTION OF POPIII OBJECTS 

In order to calculate the spatial distribution of star- 
forming halos, we employ an analytical formalism that 
tracks the correlated formation of objects. In this model, 
described in detail in Scannapieco & Barkana (2002), halos 
are associated with peaks in the smoothed linear density 
field, in the same manner as the standard Press-Schechtcr 
(1974) approach. This approach extends the standard 
method, however, using a simple approximation to con- 
struct the bivariate mass function of two perturbations of 
arbitrary mass and collapse redshift, initially separated by 
a fixed comoving distance (see also Porciani et al. 1998). 
From this function we can construct the number density 
of source halos of mass M s that form at a redshift z s at a 



comoving distance r from a recipient halo of mass M r and 
formation redshift z r : 



dn 
dMs 



{M s ,z s ,r\M r ,z r ) 



dM s dM r 



(M s ,z s ,M r ,z r ,r) 



(10) 

where -j^-{M r , z r ) is the usual Press-Schechter mass func- 
tion and d%fjM-(M s , z s , M r , z r , r) is the bivariate mass 
function that gives the product of the differential number 
densities at two points separated by an initial comoving 
distance r, at any two masses and redshifts. Note that 
this expression interpolates smoothly between all standard 
analytical limits: reducing, for example, to the standard 
halo bias expression described by Mo & White (1996) in 
the limit of equal mass halos at the same redshift, and re- 
producing the Lacey & Cole (1993) progenitor distribution 
in the limit of different-mass halos at the same position at 
different redshifts. Note also that in adopting this defini- 
tion we arc effectively working in Lagrangian space, such 
that r is the initial comoving distance between the per- 
turbations. Fortunately as enrichment is likely to be more 
closely dependent on the column depth of material be- 
tween the source and the recipient than on their physical 
separation, this natural coordinate system is more appro- 
priate for this problem than the usual Eulerian one. As a 
shorthand we define -^f(z s ,r) = Jg_(M s , z s , r\M r , z r ). 

We now wish to compute fm(M, z), the fraction of ob- 
jects of mass M forming from gas with Z < Z cr at a 
redshift z, i.e. the sites of PopIII star formation. Hav- 
ing computed the fraction of halos that are formed from 
metal-free gas, we will then be able to construct the av- 
erage number of PopIII and PopII/I outflows impacting a 
random point in space at a redshift z as 



(n 0JII (z)) = r d Zs [ 

J z J A 

(N 0JI/I (z)) = dz s 

J z J A 



dM s 



dn 



n(Zs) 

dM s 



dMAz. 
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d 2 n 
dM.dz. 
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where Vjjj = 4f r? 



3 ' iny iv±s i Zs ' z > "a — 3 
are the volumes contained in spheres of radius rui and 

Tu respectively, d ^\ n dz (M, z) is the differential Press- 
Schcchtcr number density of objects forming as a function 
of mass and redshift, and M m { n (z s ) is the minimum mass 
that can cool within a Hubble time at the specified forma- 
tion redshift z, i.e. t coo \(M Sl z s ) < tH(z s ) (see Ciardi et 
al. 2000). Finally, we include a delay between virialization 
and star formation of ts — (47rG51ol80/9 c )~ 1//2 for all ob- 
jects, which accounts for the free-fall time. This is equal 
to 2.3(1 + z) _3//2 /i _1 Gyr in our assumed cosmology. 

In constructing these integrals, we divide objects into 
two regimes, in terms of their primary cooling mechanism. 
In halos with virial temperatures above 10 4 K, atomic cool- 
ing is effective, and we adopt the fixed star formation effi- 
ciencies of fl 1 and fl 11 as described above. In halos with 
virial temperatures < 10 4 K, cooling is dependent on the 
presence of molecular hydrogen. In our previous study of 
PopII/I outflows (Scannapieco, Ferrara, & Madau 2002), 
we assumed in our fiducial model that star formation in 
these objects was completely suppressed by H 2 photodisso- 
ciation by UV radiation from the first stars. As this work 
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focuses on smaller objects, on the other hand, we allow 
star formation in such halos, but with a lower overall star 
formation efficiency. Thus for both PopIII and PopII/I 
objects with T < 10 4 K, we adopt a star formation rate 
of O.lfl 1 ' 111 . This is physically motivated by the lower ef- 
ficiency of molecular hydrogen cooling, which reduces the 
number of cooled baryons available to form stars of any 
type (Madau, Ferrara, & Rees 2001). 

With these assumptions, we can modify eq. (12) to con- 
struct (N(M r , z r )), the mean number of outflows impact- 
ing a recipient halo of mass M r and formation redshift z r . 
In this case we replace d ^ n dz (M, z) with the biased num- 
ber density of objects within the volume that is able to 
impact the recipient halo. As this number is a function of 
the distance between the source and recipient halos, Vui 
and Vn must be replaced by radial integrals. This leads 
to: 



(N(M r ,z r )) 



M min (z s ) 
12 d 2 n s 



dM s 



fin f 111 dr'r' 2 ns - r (z r') 
J in J ai i dM B dz a \ z <» ' > 



dnz 



+ { l-f ln )f^drVi^(z s y)\ll2) 
:{z s , r) is now the differential correlated num- 



where -jt-,- j- 

dM s dz y 

ber density of source objects. 

Finally, we relate this quantity to fiu(M 7 z). As our 
formalism can only account for the correlations between 
two halos, we must adopt an approximation for the frac- 
tion of objects that are affected by multiple winds. The 
simplest approach is to assume that while the winds only 
impact a fraction of the overall cosmological volume, their 
arrangement within that fraction is completely uncorre- 
cted. In this case //// and (N) share the same relation 
as they would in the uncorrelated case, 



fui(M, z) = 1 — exp[— (N(M,z))}. 



(13) 



We will adopt this Ansatz as our fiducial approach 
throughout this paper, referring to it as the "best-guess 
model" 

An alternate approximation relates these quantities by 
assuming that there is no overlap between outflowing bub- 
bles. To construct this relation, let us define /i(M, z), 
/2(M, z),... as the fraction of objects at a particular mass 
scale that are impacted by 1,2,... many outflows. In 
this case (N(M,z)} = ESi ifi(M,z) while fin(M,z) = 

i-£~i/i(^,z) = i-W«))+ai-W4 

Dropping the multiple-wind term gives the following rela- 
tion 



f in (M,z) = 



\-{N{M,z)) (N(M,z))<l 
(N(M,z))>l 



(14) 



Note that as the contribution from multiple winds is al- 
ways positive, dropping this term can only decrease the 
fraction of PopIII objects in any given model of metal dis- 
persal and enrichment. Thus this relation strictly overesti- 
mates the fraction of objects impacted by winds, and will 
be referred to as the "conservative modeV below. 

With either of these assumptions, eqs. (12) and (13) or 
eqs. (12) and (14) form a closed system, which iteratively 

1 We recall that Eg 1 ' 111 is expressed in units of lO^ergMg 1 . 



defines fjjj(M, z) and (N(M, z)). Note that as the metal- 
licity of all bubbles are well above Z CXl these equations 
are completely independent of the details of the assumed 
yields of the underlying stellar populations. 



5. THE REDSHIFT DISTRIBUTION OF POPIII OBJECTS 

Having developed an analytic formalism for PopIII for- 
mation, we now apply this method to study the evolution 
of primordial objects as a function of model parameters. 
As this evolution depends only on the spatial distribution 
of PopIII ejecta, our models can be simply classified by the 
combination of parameters entering Eg 11 , the total energy 
input into outflows per unit gas mass. To constrain this 
combination, we first consider the range of PopIII IMFs 
as compiled in Table 2. 

Surveying this table, we find that E^Af 71 <; 0.2 over 
all choices of M c and ct c , with much smaller values being 
found for the majority of models. As the star formation 
efficiency // /7 is certainly much less than 1, and is more 
likely to be on the order of 10% (Scannapieco, Ferrara, & 
Madau 2002; Barkana 2002) while f w is likely to be £ 30%, 
we can therefore place a firm upper limit of Eg 11 < 10~ 15 , 
with the most likely range being Eg 11 <; 10~ 25 (see Table 
2). 

In the left panels of Figure 2 we plot the fraction of 
PopIII objects, fiu(z), in the best-guess model 1 with 
Eg 11 = 10~ 15 as a function of mass, as well as for weaker 
feedback cases in which E[ u = 10~ 2 - 5 , 10~ 3 - 5 , and 10" 4 5 . 

In each model we have taken fjM 11 = 4.8 x 1O _4 M ^ 1 , 
consistent with models SNII-B and SNII-C with a fixed 
star formation efficiency of = 0.1, and have calculated 
eqs. (10) - (13) for objects at 25 different mass scales, 
spaced in equal logarithmic intervals from 1.0 x 10 6 M Q 
to 1.0 x 10 n M Q . 

In all cases, the large masses are the most clustered and 
thus feedback affects them more strongly, resulting in a 
quicker transition from primordial to pre-enriched star for- 
mation. In the Eg 11 = 10~ 15 model, in fact, only objects 
with masses ~ lx 1O 7 M are able to form PopIII stars, 
while at higher mass scales the strong feedback of this 
model suppresses primordial star formation almost imme- 
diately. At smaller mass scales, on the other hand, the 
minimum mass cutoff occurs at very early times, resulting 
in the sharp drop seen in the 1O 7 M curve. Thus it is 
only in a narrow range of masses and redshifts that PopIII 
objects are able to be formed in the presence of strong 
feedback. 

Reducing the energy input helps to ease these restric- 
tions, as can be seen in the Eg 11 = 10~ 2 ' 5 and 10~ 3 ' 5 
models. In the latter case there is a significant population 
of objects with masses of a few times 10 8 M Q , while the 
~4x 10 7 M Q scale objects continue to be formed down 
to their 10 4 K limit at z ~ 5. Even in this model, how- 
ever, the majority of PopIII objects tend to be in the 
10 6 ' 5 M Q — 1O 7 '°M0 mass range, just large enough to cool 
within a Hubble time, but small enough that they are not 
clustered near areas of previous star formation. 

At even smaller feedback values, the mass and redshift 
ranges of PopIII formation widen to the point that they 
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conflict with observations. Thus in the E 1 ^ 1 = 10~ 4 ' 5 
model, metal dispersal is so inefficient that the universe 
is unable to break out of the star-formation conundrum. 
Pair production SNe are too infrequent to spread material 
over a large fraction of the IGM, and PopIII star formation 
continues indefinitely. 

The same trends are seen in the conservative models, 
plotted in the right panels of this figure. While feed- 
back is more severe in this case, reducing E^ 11 again re- 
sults in a systematic widening of the mass and redshift 
range of PopIII formation. As in the previous case, at 
the lowest £^ n value, the cosmological enrichment be- 
come extremely inefficient, and universe become ensnared 
in the star-formation conundrum. As this model is con- 
structed to systematically underestimate the fraction of 
primordial objects, this result can be taken as a conser- 
vative lower bound on feedback. Thus we can restrict 
our attention to feedback models with energy inputs the 
range 10~ 15 > E^ 11 > 10 -4 ' 5 , with a most likely range of 

io- 2 - 5 £ e{ u % 10- 4 . 

In the left panels of Figure 3 we plot the overall star 
formation rate per comoving Mpc 3 in PopIII and PopII/I 
objects, for two choices of fl 11 ■ We compute these values 
simply as 



fin 



-1 cIII 



f-lC 



and we show our results 
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dtdM 



{O.1 + O.90[M - M 10 4(i)]} , (16) 



where 9 is the Heaviside step function, which is used to 
account for the O.lfJ. 11 star formation efficiency assumed 
objects with masses smaller than M W i(t), the mass cor- 
responding to a virial temperature of 10 4 if at the time 
t. The overall downturn at z ~ 4 in these plots is due 
to fact that we are only considering burst-mode star for- 
mation, combined with our choice of an upper mass limit 
of 1O 12 M , which crudely models the lack of larger star- 
bursts due to the long cooling times in these objects. Be- 
cause quiescent- mode star formation in Pop II/I objects 
becomes important below this redshift, we restrict the 
analysis given in Figs. 3 to 8 to rcdshifts above 4. Note that 
in all cases PopIII star formation in objects with masses 
^ 1O 1O M is negligible, and we impose fm(M,t) = for 
all objects with masses greater than 10 11 M Q , as we do not 
calculate these values directly. 

While the peak of PopIII star formation occurs at z <~ 10 
in all models, such stars continue to contribute apprecia- 
bly to the SFR density at much lower redshifts. This is 
true even though the mean IGM metallicity has moved 
well past the critical transition metallicity. In the right 
panels of Figure 3, we show the average values as com- 
puted in our models. In this case, rather than fix 
we must now estimate an overall number of SN 77 per unit 
mass in PopIII objects for each model. To do this, we 
assume as in Figure 1, that Y mct /M Q ~ 2£kin/10 51 ergs, 
and that half the metals are ejected into the IGM. In this 



case Y met N^fi' _ n , 
for models in which f w = 0.3 and /„, = 0.1. Note that 
the metallicities are lower in the f^ 11 = 0.3 models than 
in the f^ 11 = 0.1 models, as increasing f w while leaving 
Eg 11 fixed corresponds to a decrease in the overall number 
of SN 77 in each object. In all cases, however, the aver- 
age IGM metallicity exceeds Z cr at very early redshifts of 
12 — 20, consistent with previous estimates by Schneider 
et al. (2002) and Mackey, Bromm, & Hernquist (2002). 

Unlike these previous investigations, however, our cor- 
related structure formation model indicates that it is not 
until much lower redshifts that outflows are able to spread 
enriched gas to all regions in which galaxy formation is 
taking place. This inefficiency of ejection means that ap- 
preciable levels of PopIII star formation arc likely to con- 
tinue well into the observable redshift range of z <J 6, con- 
tributing to as much as 10% of the star formation rate at 
z = 6, for our most favorable choice of parameters. Thus 
it is likely that the direct observation of PopIII objects is 
well within the capabilities of current instruments, and in 
fact such objects may have already been seen in ongoing 
surveys of high-redshift objects. 

6. THE DIRECT DETECTION OF POPIII STARS 
6.1. Properties of PopIII Objects 

Our findings in §5 have important implications for the 
development of efficient strategies for the detection of 
PopIII stars in primeval galaxies. This is because at any 
given redshift, a fraction of the observed objects have a 
metal content that is low enough to allow the preferential 
formation of PopIII stars. As metal-free stars are power- 
ful Lya line emitters (Tumlinson, Giroux, & Shull 2001; 
Schaerer 2002; Venkatesan et al. 2003), it is natural to use 
this indicator as a first step in any search for primordial 
objects. This is even more promising as surveys aimed at 
finding young, high redshift systems have already discov- 
ered a considerable number of such emitters (Dey et al. 
1998; Weyman et al. 1998; Hu et al. 1999, 2002; Ajiki et 
al. 2002; Dawson et al. 2002; Frye, Broadhurst, & Bem'tez, 
2002; Lehnert & Bremer 2002; Rhoads et al. 2002; Fujita 
et al. 2003; Kodaira et al. 2003) 

In principle, the calculation of the Lya luminosity, L a , 
from a stellar cluster is very simple, being directly related 
to the corresponding hydrogen ionizing photon rate, Q(H), 

by 

L a = c L (l- f csc )Q(H), (17) 

where cl = 1-04 x 10~ n erg and / esc is the escape frac- 
tion of ionizing photons from the galaxy. As this escape 
fraction is relatively uncertain, we adopt here an educated 
guess of / osc = 0.2, which is based on a compilation of both 
theoretical and observational results (see Ciardi, Bianchi, 
& Ferrara 2002 for a detailed discussion and references 
therein). Finally, Q(H) is computed from evolutionary 
stellar models. In this case, we assume that in all objects, 
the stars giving rise to the observed Lya emission form in 
a burst, which is an increasingly reasonable assumption for 
higher redshift /lower mass systems. We take this burst to 
be coeval with the galaxy formation redshift for all objects. 
For PopIII stars, we rely on the recent work of Schaerer 
(2002) who computed the time evolution ofQ(H) in such a 
burst for several assumed IMFs. We consider two models, 
assuming that that these stars are distributed according 
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to a Salpeter IMF with a low-mass cutoff equal to either 
M, = 1M ("normal") or M t = 5OM ("top-heavy"), and 
an upper mass cutoff held fixed at M u — 500M Q . For 
PopII/I stars, on the other hand, we make use of the stel- 
lar evolution code STARBURST99 (Leitherer et al. 1999), 
assuming in all cases that M; = 1M Q and M u = 120Mq 
with a metallicity equal to 1/20 Zq. 

Fig. 4 shows the evolution of the ionizing photon rate as 
a function of the time elapsed since the burst in these mod- 
els. As expected, PopIII stars are characterized by shorter 
lifetimes and much higher luminosities, principally due to 
their higher surface temperatures. Also note the sudden 
drop of Q(H) at about 40 Myr, corresponding to the life- 
time of stars more massive than w 8M Q . For comparison 
we have also plotted Q(H) for a 6(M+ = 260 M Q )-function 
IMF, which provides an extreme example of a bright, short 
Q(H). These curves were directly applied to calculate the 
Lya emissivity of PopIII and PopII/I objects as a function 
of their stellar mass and redshift according to eq. (17) and 
the formalism introduced in §5. It is important to note 
that in adopting this approach, we are implicitly ignoring 
the destruction of Lya photons both by dust and by the 
random motions of the medium through which they prop- 
agate (Neufeld 1991; Haiman & Spaans 1999). To calcu- 
late these effects, one must consider a number of physical 
properties that are difficult to model, including metallicity, 
dust content, gas random motion and column-density dis- 
tribution. Although this has been attempted by Haiman 
& Spaans (1999), we prefer here to avoid these intrica- 
cies and consider our values as upper limits to the actual 
Lya emission. 

Having calculated the Lya emission, we can then derive 
the probability that a given high-rcdshift Lya detection 
is due to a cluster of PopIII stars. Here we allow for the 
time evolution of the Lya line luminosity from the burst 
down to the redshift at which the object is observed ac- 
cording to the curves in Fig. 4. Moreover, we assumed, 
as in eqs. (16) that fj 1 = fj 11 = 0.1 for all objects with 
masses greater than M 10 4(t) and f^ 1 = fl 11 = 0.01 for 
masses between M m i n (i) and M W i(t). Finally, we assume 
a top-heavy IMF for PopIII stars, which greatly shortens 
the average lifetimes of PopIII systems and boosts their 
luminosities by almost an order of magnitude. The result- 
ing probabilities are displayed in Fig. 5. In this figure, the 
isocontours in the Lya luminosity-redshift plane indicate 
the probability to find PopIII objects (i.e. galaxies hosting 
PopIII stars) in a given sample of Lya emitters for various 
feedback efficiencies, parameterized by the value of £^ n . 
For reference, we also include the available data points (as 
described in the caption) and a line corresponding to a 
typical detection flux threshold of 1.5 x 10~ ergs cm~ 2 
s _1 . For the Fujita et al. (2003) data, we have included a 
random scatter of Az — ±0.22 to simulate observational 
uncertainties. 

In this figure, we see that PopIII objects populate a 
well-defined region of the L Q -redshift plane, whose extent 
is governed by the feedback strength. Note that the lower 
boundary is practically unaffected by changes in Sg 11 , as 
most PopIII objects are in a limited mass range such that 
they are large enough to cool efficiently, but small enough 
that they are not clustered near areas of previous star for- 
mation. At lower £^ n values, the non-zero probability 



area widens considerably: in this case, a smaller volume 
of the universe is polluted and PopIII star formation con- 
tinues at lower redshifts and in the higher mass, more lumi- 
nous objects that form later in hierarchical models. Above 
the typical flux threshold, Lya emitters are potentially de- 
tectable at all redshifts beyond 5. Furthermore, the frac- 
tion of PopIII objects increases with redshift, independent 
of the assumed feedback strength. For the fiducial case 
£g 77 = 10~ 3 , for example, the fraction is only a few per- 
cent at z = 4 but increases to approximately 15% by z = 6. 
We then conclude that the Lya emission from already ob- 
served high-z sources can indeed be due to PopIII objects, 
if such stars were biased to high masses. Hence collecting 
large data samples to increase the statistical leverage may 
be crucial for detecting the elusive first stars. 

The situation is markedly different, however, if the first 
stars formed according to a normal IMF, as shown Fig. 6. 
In this case, the chances to detect PopIII stars at low red- 
shifts rapidly decreases as their flux drops below the sen- 
sitivity threshold. This is particularly evident when the 
feedback is strong and cosmic metal enrichment proceeds 
efficiently. In this case, detections are far more difficult, 
although an observed absence of galaxies hosting PopIII 
stars could be used to set lower-limits on the efficiency of 
feedback processes in the early universe. 

Our simple burst model of high-redshift sources can also 
be compared to the observed number counts of Lya emit- 
ters per unit redshift and square degree. In Fig. 7, we plot 
these number counts for both for PopII/I objects, which 
are largely independent of E 1 ^ 1 , and PopIII objects whose 
numbers depend sensitively on this value. These curves 
are computed in both the top-heavy and normal IMF cases 
and arc contrasted with the data by Hu et al. (1998). In all 
cases PopII/I objects dominate the number counts at red- 
shifts <; 6, However, at higher redshift the PopIII objects 
outnumber "normal" galaxies for a heavy IMF, provided 
the feedback is mild. In all cases the detection probability 
increases rapidly with redshift, as we saw above. 

Finally, we show in Fig. 8 the cumulative probability 
distribution function (PDF) of observed Lya equivalent 
widths (EW) characterizing PopIII objects that fall above 
the assumed detection threshold of 1.5 x 10 ergs cm~ 2 
s _1 . Here we consider a top-heavy IMF and two models 
assuming a l-500M o Salpeter IMF with \og(S^ n ) = -3.5 
and log(£ 7// ) = -4.0. Note that the top-heavy IMF case 
is independent of both redshift and £^ n because the Lya 
luminosity is roughly constant over the short lifetime of 
such massive stars, and thus the objects remain detectable 
over the full evolution of equivalent widths. For all models, 
the equivalent width PDF docs not depend dramatically 
on the assumed IMF or the redshift, and PopIII stars are 
able to populate the large EW (log EWl Vc , > 3) region of 
the graph, although such objects are observed (eg. Malho- 
tra & Rhoads 2002). This region is impossible to fill with 
PopII/I stars, as these stars tend to inflate the low EWs 
tail of the PDF. Hence, large EWs are a key indicator of 
the presence of primordial stars. 

6.2. Strategies for Detection and Identification 

The properties identified above can be combined to yield 
a strategy to detect the cosmic objects hosting the first 
stars. Obviously, this would be a spectacular discovery, as 
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local searches for truly metal-free stars have been painfully 
unsuccessful (for a review see Beers 2000). We propose 
here to search for these stars in high-z Lya emitters, where 
chances might be considerably higher. 

The individuated strategy would be as follows: 

• Collect a large sample of Lya emitters, preferably 
weighted to the highest possible redshifts. 

• Among the objects in this sample, make a priority 
list ranking the possible PopIII candidates on the ba- 
sis of their luminosity and redshift according to the 
probability distributions shown in Figs 5 and 6. This 
can be done modulo the choice of an IMF for PopIII 
stars and a value of the feedback strength, but it is 
not overly dependent on this choice. In general the 
highcr-rcdshift objects, close to the detection limit 
are favored in most models. 

• Further select the candidates on the basis of their 
Lya EW, sources showing the largest values receiv- 
ing the highest ranking. 

• Return to the most probable candidates and perform 
extremely deep, long exposure spectroscopy. 

The goal is then to find a clear-cut feature that uniquely 
identifies a given high-rcdshift stellar cluster as virtually 
metal-free. Here they key indicator may have been pointed 
out by Tumlinson, Giroux, & Shull (2001), who noted that 
Population III stars produce large He III regions, which 
can emit detectable He II recombination emission. Re- 
combination lines of He II at A1640 A(n = 3 -» 2), A3203 
A(n = 5^ 3), and A4686 A(n = 4^3) are particularly 
attractive for this purpose because they suffer minimal ef- 
fects of scattering by gas and decreasing attenuation by 
intervening dust. Fig. 1 of Tumlinson, Giroux & Shull 
(2001) shows the flux of the He II 1640 line as a function 
of source redshift. For star formation rates of 20 M yr _1 
and 5 M yr _1 , the A1640 flux is detectable out to z w 5 
at the sensitivity level of current Lya emission line sur- 
veys, although the flux for A4686 is 7.1 times lower. These 
lines arc uniquely produced with this intensity by massive 
PopIII stars. 

As an example, let us consider the particular case of 
the highly (w 30 times) magnified Lya emitter detected 
by Ellis et al. (2001), located at redshift z = 5.576. This 
system has as an estimated stellar mass of 10 6 M Q and a 
current star formation rate derived from the 2 x 10~ 18 ergs 
cm~ 2 s _1 Lya line of 0.5 M yr _1 , implying an approx- 
imate age of 2 Myr. No stellar continuum is detected to 
an upper limit of 3 x 10~ 20 erg cm~ 2 A -1 . For this ob- 
ject, the He II 1640 emission is redshifted roughly into the 
Z band (1.078 /zm) where the difficulties of ground-based 
NIR spectroscopy make detection challenging, but not im- 
possible. In this case, the object falls between the OH sky 
lines, but it is likely to be very faint. The relative flux ra- 
tio, T of the Lya to He II A1640 lines is equal to 17.8/^j, 
where 0.4 < f evo i < 2.0 is a parameter which accounts for 
the time evolution of the stellar ionizing continuum radia- 
tion. Thus the He II A1640 line is about 9-45 times fainter 
than Lya , corresponding to 4-20 xlO -20 ergs cm' 2 s" 1 , 
a difficult experiment, but within in the realm of possibili- 
ties for some existing/planned instruments. As many high 
redshift sources are over an order of magnitude brighter, 



however, the discovery space for metal-free galaxies is large 
and deserves continued observational emphasis. 

• 7. POPIII SIGNATURES IN THE INTRACLUSTER MEDIUM 

Our PopIII models can also be directly applied to asses 
the impact of the first stars on galaxy clusters, the largest 
virialized structures in the universe (see Rosati, Borgani, 
& Norman 2002 for a recent review) . In hierarchical mod- 
els of structure formation, clusters are predicted to form 
from the gravitational collapse of regions of several Mpc, 
associated with rare high peaks in the primordial density 
field. In the process of virialization, a hot diffuse gas is 
formed, the intracluster medium (ICM), which represents 
a substantial fraction of the cluster mass, ~ 15%. The 
ICM permeates the cluster gravitational potential well and 
emits in the X-ray band via thermal bremsstrahlung. 

Rich clusters of galaxies can be considered to be "closed 
boxes" , isolated systems that reflect the overall history of 
structure formation with negligible interaction with the 
surroundings. In particular, the ICM should maintain 
clear imprints of the thermal and chemical evolution of the 
baryons and is therefore a promising location to look for 
signatures of PopIII stars (see also Loewenstein 2001). In 
the following sections we briefly review some of the avail- 
able observations of the ICM obtained by the old (ASCA, 
Bcppo-Sax) and current generations of X-ray satellites 
(Chandra, Newton-XMM). Then we present the predic- 
tions of the model that we have developed in the previous 
sections and finally discuss our results and compare them 
with previous analyses. 

7.1. ICM Abundances and Heating: Observational 
Constraints 

The hot ICM has been extensively studied via X-ray 
imaging and spectroscopy, and is observed to be enriched 
to a significant fraction of solar metallicity (see e.g. Ren- 
zini 1997; White 2000 and references therein). Although 
stars represent only about 10% of the total baryonic mass 
of clusters of galaxies (Loewenstein 2000), the favored 
mechanism for enriching the ICM is supernovae-driven 
winds from elliptical member galaxies. Various galactic 
wind models have been proposed (see e.g. Fusco-Femiano 
& Matteucci 2002; Pipino et al. 2002) but whether these 
winds are dominated by the ejecta of SNII or SNIa is 
still vigorously debated (Tsujimoto et al. 1995; Gibson, 
Loewenstein, & Mushotzky 1997; Finogucnov, David, & 
Ponman 2000; Gastaldello & Molendi 2002). SNIa prod- 
ucts are iron-rich whereas those of SNII are rich in a- 
elements, such as Si and O; therefore it is crucial to have 
accurate separate measurements of each of these contribu- 
tions. 

In the following we consider only average elemental 
abundances and neglect the detailed spatially-resolved 
abundance measurements that can be obtained with Chan- 
dra and XMM. Our reasoning is that the present study is 
best suited to make average predictions rather than appli- 
cations to specific systems. We consider elemental abun- 
dances that have been averaged over ensembles of clusters 
with similar properties. Abundance measurements are al- 
ways normalized to solar, with solar photospheric values 
from Anders & Grevesse (1989), and refer only to the outer 
regions of clusters to avoid contamination by central gra- 
dients. 
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A compilation of the available data is shown in Fig. 9. 
Triangles indicate ASCA observations and are taken from 
Table 2 of Fukazawa et al. (1998). Dots indicate the results 
of a recent analysis of ASCA archival data by Baumgart- 
ner et al. (2002). Thanks to the high spectral resolution 
and better sensitivity, XMM-Newton has obtained the first 
abundance measurements of oxygen. Recently, Tamura et 
al. (2002) have reported new results of abundance mea- 
surements and have consistently analyzed all the available 
data. The sample consists of groups and poor clusters, 
with temperatures in the range [1.2 - 4] keV. The abun- 
dance ratios in the outer regions of the observed clusters 
are shown in Fig. 9 as solid regions. The results for Si 
and S are consistent with other measurements but there is 
a significant change in the O/Fe abundance between the 
center and the outer region, varying from 0.5 up to 0.8 
solar. 

We show in Fig. 10 the relative abundances averaged 
over three temperature bins, [2 - 4] keV, [4 - 8] keV, [8 
- 12] keV 2 . Unfortunately, measurements of the (O/Fe) 
abundance are restricted only to the first bin. However, 
the relative abundances of the other elements show some 
trends: going from the poorer to the richer systems, Fe 
decreases whereas (Si/Fe) and (S/Fe) increase. This is 
consistent with previous findings (Fukazawa et al. 1998; 
Finogucnov et al. 2000) and might reveal a larger role of 
SNIa in the enrichment of groups and poorer systems com- 
pared with that of rich clusters. Supersolar (Si/Fe) and 
(S/Fe) abundance ratios are observed in richer clusters. 

Besides elemental abundances, the observed X-ray prop- 
erties of clusters provide strong constraints on the thermal 
history of the ICM. While the optical properties of clus- 
ters can be understood in the context of self-similar mod- 
els (Kaiser 1986), which assume pure gravitational heat- 
ing, the slope of the ICM X-ray luminosity-temperature 
(L x — T) relation hints a more complicated history. If 
the ICM was heated only gravitationally, then the X-ray 
luminosity should be proportional to T 2 , yet observations 
indicate that this luminosity instead goes as T 7 ' 2 . At lower 
temperatures this relation become even steeper, with clus- 
ters of widely different luminosities all having tempera- 
tures ~ 1 keV (Cavaliere, Menci, & Tozzi 1999). 

The accepted explanation for this unusual behavior re- 
lies on the presence of large-scale energy input from astro- 
physical sources (Kaiser 1991) prior to the gravitational 
collapse of the cluster. This energy increases the ICM en- 
tropy, placing it on a higher adiabat and preventing it from 
reaching a high central density during collapse, which de- 
creases its X-ray luminosity (e.g. Tozzi & Norman 2001). 
For a fixed degree of heating per gas particle, this effect is 
more prominent for poorer clusters, whose virial temper- 
atures are comparable to this extra contribution. As a re- 
sult, a T 7 / 2 relation is established in hot systems and bro- 
ken for colder systems, as observed. Both semi-analytical 
studies (e.g. Cavaliere, Menci, & Tozzi 1998; Wu, Fabian, 
& Nulsen 2000; Tozzi, Scharf, & Norman 2001) and numer- 
ical simulations (e.g. Brighenti & Mathews 2001; Bialck, 
Evrard, & Mohr 2001; Borgani et al. 2001) indicate that 
~ 1 keV per gas particle of extra energy is required, with 
a corresponding entropy floor of 55 — 110 h~ 1/>3 keV cm 2 
(see however, Voigt & Bryan 2001). 



7.2. Predictions for PopIII stars in clusters 

Using the distribution of PopIII objects calculated in 
§3 we can directly relate E^ 11 to the levels of enrich- 
ment and heating observed in clusters. In order to com- 
pare these quantities in a model-independent way, we 
work with Fm , the fraction of cluster gas that cools into 
PopIII objects. While this quantity is dependent only on 
the quantity £g 7/ , it can be easily used to construct the 

overall fraction of PopIII stars (fj 11 Fm), energy input 
( £ ni Nllf ni Fiii ^ and metal yidds (Y^N^fl H F ni ) 

in any given model of PopIII star formation. 

Accounting for the 0.1/ 7// efficiency assumed for objects 
with T v j r < 10 4 K, we calculate this quantity as 

Fm = ST dz *rt ImIm dM s M sflII (M, z s ) 



dn s 



' dz s dM. 



{0.1 + 0.98[M S - M 10 4z s )}} (18) 



In this case the relevant recipient halo is the cluster, and 
all sources are located within the region of space that later 
collapses into the cluster itself. As ( z c r ) coincides 

almost exactly with fj 1 ^ (z r ,0) if the distance between 
the two peaks is less than their collapse radii, we are free to 
take the number density at zero separation (Scannapieco 
& Barkana 2002), which is equivalent to the usual progen- 
itor distribution as first described in Lacey & Cole (1993). 
While these integrals break down for the smallest objects 
and redshifts, the rising value of M m i n , coupled with the 
sharp drop in fm(M, z) at low redshifts assures that the 
contribution from these values is minimal. Nevertheless we 
impose an overall minimum source redshift of 4 in comput- 
ing this integral, rather than include contributions down 
to the final collapse redshift z c . 

In Figure 11 we plot Fm(£^ n , z c ) assuming both the 
best-guess (solid) and conservative (dashed) models. In 
both cases we consider clusters with collapse redshifts of 
2.0, 1.5, 1.0, and 0.5. Note that these values are completely 
independent of the mass of the final cluster, which is typi- 
cally many orders of magnitude larger than the masses of 
PopIII objects. From the point of view of PopIII objects, 
the regions of space that form clusters are equivalent to 
closed universes in which the overall matter density has 
been increased such that collapse occurs at z c . 

A clear implication of this plot is that at most 10% of 
the gas in clusters can be cooled into PopIII objects. This 
is true even in the most optimistic case in which £^ n is 

taken to be 1.0~ 4 ' 5 M Q -1 and the best-guess model, eq. 
(13), is used to relate //// and (iV). In the more likely 
range of values 1.0~ 40 <^ E^ 11 <^ 1.0~ 25 , Fjjj varies from 
about 1% to 5% of the cluster gas. When combined with 
the detailed PopIII modeling discussed in §2.1, these lim- 
its place strong constraints on the PopIII contribution to 
ICM heating and enrichment. 



7.3. ICM Heating 

Our Fm values can be simply related to the constraints 
on cluster preheating without any further assumptions 
about the underlying distribution of stars. In order to 

2 For each temperature bin, we have averaged over data points obtained by different groups as shown in Fig. 9, and we have associated only 
the statistical error with this average, neglecting any systematic uncertainties. 
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place a firm upper limit on this contribution, we assume 
that while only a fraction f w of the kinetic energy from 
SN 77 goes into powering outflows and dispersing metals, 
the rest of the energy also somehow makes its way into the 
ICM. In this case a total energy input of 

£i"N^fl n F nl crgs/M = aiOf-iE*" Fin kev/particle 

(19) 

is injected into this gas. Finally we take a conservative 
value oi f w = 0.1 for the overall wind efficiency, such 
that the energy input into the ICM is large relative to 
the outflow energies. In the lower panel of Figure 11, we 
plot this estimate for both the best-guess and conserva- 
tive models. Focusing our attention on the upper curves, 
we see that while increasing £ 777 results in a higher over- 
all level of heating, the energy input per baryon is always 
Ticm < 0.2 keV/particle, even for the maximum reason- 
able £ 777 value of 10 -15 . Thus SN 77 heating falls short of 
the required level of ~ IkeV per gas particle even in the 
most optimistic case. 

Furthermore, most of this heating occurs at high red- 
shifts, at which the mean cosmological density is higher 
and thus the overall entropy corresponding to a given tem- 
perature is smaller. Taking a rough value of the redshift for 
this energy injection of z = 8 (see Fig. 3) and ignoring the 
additional gas density due to the fact that clusters form 
from overdense regions gives an overall entropy level of 
2/3 - 150/i- 1/3 7i CM kev cm for our assumed cos- 
mology. Even in the most optimistic heating model, this 
is barely consistent with the lowest allowed entropy values 
of 30/i~ 1//3 kcV cm 2 . Thus we conclude that the energy 
input from PopIII objects is insufficient to preheat clus- 
ters to their observed levels. It is important to emphasize 
that this result is based simply on energetic arguments, 
and does not depend on any assumptions about the un- 
derlying PopIII star formation rate, initial mass function, 
or properties of SN 77 . 

7.4. ICM Metal Enrichment 

In order to estimate the observed (Fe/H), (Si/Fe), 
(S/Fe) and (O/Fe) abundances from our models we com- 
pute the total mass of a given element i in the ICM as 



M l ^Af n M,->Y/ I +Y l 



la 



Kfll 

y.TT_ n 

1 M 11 



(20) 



where Y/ 1 , Y/ a and Y? 1 are the IMF- averaged yields of 
the element i (in solar masses) for SNII, SNIa and SN 77 as 
given in Table 3, N 11 and A/" 77 are the number of SNII 
and SN 77 per stellar mass formed as given in Tables 1 
and 2 and the relative number of SNIa and SNII is left 
as a free parameter, which is later constrained by obser- 
vations. Finally, M* is the total mass of PopII/I stars 
formed assuming a Salpeter IMF and 1Z = M/ 77 /M* is 
the total mass of PopIII stars formed per solar mass of 
PopII/I stars, assuming a PopIII IMF as described in sec- 
tion 2.2. The abundance ratio of a given element i relative 
to iron can be expressed as 



(i/Fe) 



Y/ 1 + Y/ a (Af Ia /Af H ) + YP{N^/N u )n 
y fo + Y^{N Ia IN H ) + Y^(N^/N H )1Z : 



(21) 



and therefore depends only on the relative rate of SNIa 
over SNII and on the parameter 7Z. Conversely, the abso- 
lute Fe abundance is computed as follows, 



(Fe/H) 



Mi 



Fe 



Mi 



Fc 



0.1M Fe 



O.lMpc 



M H 0.76M gas 0.76M* 0.76(1 - i?)M* 

(22) 

where we have assumed that (i) 76% of the gas mass is in 
hydrogen (ii) the present-day mass in stars and remnants 
M^o is 10% of the gas content (in) we have corrected for 
the fraction R of the stellar mass of PopII/I stars which is 
returned to the gas. Assuming a Salpeter IMF and rem- 
nant masses from Ferreras and Silk (2000), we find that 
R = 0.276. Thus, the absolute iron abundance is, 



(Fe/H) = 



0.LV 



ii 



0.76(1 - R) 



AT 



A/" 77 

W 7 



K 



(23) 

Our modeling of Fin gives us a self-consistent estimate 
for the parameter 1Z. In particular, the parameter 1Z can 
be related to Fni as, 



1Z 



M 



/// 



fl n M™ 



I. 



in 



M* 



(l-_R)M* i0 0.1(1 -R) 



F, 



in- 



(24) 



For a given PopIII IMF and efficiency / 777 , we can esti- 
mate £g 77 , taking a conservative value of f w = 0.1 (see 
the last column of Table 2), and find the corresponding 
values for Fm and 1Z. The results are presented in the 
top panels of Fig. 12, where we show 1Z as a function of 
£ 777 for two possible cluster formation redshifts (0.5 and 
2) for two different choices of the PopIII star formation 
efficiency (0.1 and 0.5). Each set of horizontal lines indi- 
cates the predicted £ 777 for SN 77 -B (dashed lines), SN 77 - 
C (solid lines) and SN 77 -D (dot-dashed lines) assuming 
(Tp ln and cr™ ax . The intersection points mark the appro- 
priate 1Z for each model. It can be seen from the figure 
that 1Z is typically between 1% and 2%, and never exceeds 
5%. This result is very solid as 1Z is largely independent 
of the assumed cluster formation redshift and of the de- 
tailed PopIII IMF. Given the range of values for 1Z in each 
SN 77 model, we can estimate the corresponding Fe abun- 
dance normalized to solar using eq. (23) and assuming that 
only PopIII stars contribute to the ICM metal enrichment. 
The results are shown in the bottom panels of Fig. 12 us- 
ing the same line coding as the top panels. The predicted 
Fe abundance in model SN 77 -C is degenerate in cr™ m and 
cr™ ax and therefore appears as a solid thick line instead of 
as an extended region. This figure shows that the typi- 
cal [Fe/H] abundance that can be contributed by PopIII 
stars is at best ~ 0.07. Since the observed iron abundance 
ranges between 0.2 and 0.4 solar, this means that PopIII 
stars can contribute no more than 17%-35% of the iron 
observed in clusters. 

The predicted elemental abundance ratios are shown in 
Figs. 13 and 14 as a function of the relative number of 
SNIa to SNII. The shaded regions in each panel indicate 
the range of observed values for (Fe/H), (Si/Fe), (S/Fe) 
and (O/Fe) averaged over the same three temperature bins 
as in Fig. 10. In Fig. 13 we show the predicted abundances 
if only SNIa and SNII enrich the ICM, i.e. if 1Z = 0. The 
different lines correspond to different progenitor models for 
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SNII (see Table 1). It is clear that none of these models 
can simultaneously account for the observed abundances, 
as previously recognized by Loewenstein (2001). This is 
true even if one assumes different relative rates of SNIa and 
SNII depending on the cluster richness. Recently Tamura 
et al. (2002) claimed that XMM observations of poor clus- 
ters were consistent with enrichment by SNIa and SNII 
with relative rates in the range 0.2-0.7. We find, however, 
that when averaged with ASCA observations at compara- 
ble richness, the agreement is no longer valid as it over- 
predicts (S/Fe). 

Hereafter, we take as our fiducial model SNII-C and con- 
sider whether PopIII stars with their peculiar nucleosyn- 
thetic yields can improve the agreement with the observed 
data. The results are shown in Fig. 14 where we add a 
contribution from PopIII stars formed according to model 
SN 77 -B with <r™ m . We restrict our analysis to this model 
as it provides the best fit to the observed abundances. In 
each panel, abundances are computed assuming the lim- 
iting values for 1Z consistent with the model (see Fig. 12) 
and are compared with the prediction of model SNII-C 
(solid line). 

Because of the relatively small contribution of PopIII 
stars to the total iron abundance (see Fig. 12), N Ia /N 11 
values from 0.3—1 are required to match the observed data, 
with higher values corresponding to poorer clusters. How- 
ever, pre-enrichment by PopIII stars enhances the Si/Fe 
ratio and the same relative rate of Type la to Type II SNc 
can simultaneously account for (Fe/H) and (Si/Fe). Since 
S and Si trace each other in SN 77 ejecta, the increase in the 
Si/Fe ratio is reflected in a similar increase in S/Fe which 
tends to be systematically overpredicted with respect to 
the observed data. Finally, the O/Fc abundance is always 
underestimated for the relative Type la over Type II SN 
rates required to match Fe and Si/Fe data being 0.4 so- 
lar for model SN 77 -B. This is roughly consistent with the 
value observed by XMM in the inner cluster regions. 

To summarize, PopIII stars can contribute no more than 
20 % of the iron observed in galaxy clusters but if they 
form with characteristic masses of about 200 — 260M Q , 
their peculiar elemental yields help to reconcile the ob- 
served Fe and Si/Fe abundances. However, they tend to 
overpredict S/Fe and can account only for the O/Fe ratio 
in the inner regions of poor clusters. Additionally, their 
energy input is insufficient to heat the ICM to the observed 
levels. 

Observations of elemental abundances will increase in 
the near future, thanks to the enhanced sensitivity of the 
new generation of X-ray satellites. The presently available 
data appear to be hard to reconcile with simple enrich- 
ment by SNIa and SNII ejecta and are only marginally 
improved by the addition of a phase of pre-enrichment by 
very massive PopIII stars. The observed iron abundance 
requires, for poor clusters, very high relative rates of Type 
la over Type II SNe, at least if their progenitors are formed 
according to a Salpeter IMF. 

8. CONCLUSIONS 

Although the formation of metal-free stars marked the 
dawn of the modern universe, their properties remain a 
cosmological mystery. As the cooling rates of collapsing 
clouds are drastically changed by even trace amounts of 
metals, fragmentation in primordial regions is completely 



different than those in molecular clouds today. Similarly, 
the processes that regulate protostellar accretion rates in 
PopII/I stars failed to operate in primordial gas, radically 
altering evolution. Thus theoretical PopIII star-formation 
studies, while in many ways more straightforward than 
studies of later-forming stars, can not be directly com- 
pared with any observational constraints. And although 
reduced cooling and lack of regulation of accretion both 
point to massive star formation, the characteristic mass of 
such stars is uncertain to over an order of magnitude. 

Observations of PopIII stars have remained equally elu- 
sive. Not a single metal-free star has been seen in our 
Galaxy, lending weight to the top-heavy hypothesis while 
at the same time forcing observers to adopt alternative 
methods. These include relating the elemental abundance 
patters observed in metal-poor halos stars with the prop- 
erties of the (unobserved) stars that enriched them, an ap- 
proach that has uncovered unusual abundance ratios that 
may point towards SN 77 enrichment. A second method 
is to examine the metal content of the high-redshift IGM 
itself, although such studies are limited by the small num- 
ber of observable elements and uncertain ionization cor- 
rections. 

In this work, we have explored two alternative meth- 
ods for studying PopIII stars: (1) by direct detections of 
high-redshift sources, which are accessible to ongoing sur- 
veys of Lya emitters; (2) by studying the impact of PopIII 
stars on the intraclustcr medium, which has been well ob- 
served in the X-ray. As the properties of the first stars 
are extremely uncertain, our investigation has relied on 
analytical models that allow us to make the most general 
statements possible based on a minimum number of as- 
sumptions. 

Using a simple 1-D outflow model, we have found that 
for the full range of parameters considered, the mean 
metallicity of outflows from PopIII objects is well above 
the critical threshold that marks the formation of normal 
stars. This conclusion is based only on three assumptions: 
(1) that Ymot/M© <~ 2£kin/10 51 ergs for SN 77 as computed 
by Heger and Woolsley (2002); (2) that the critical transi- 
tion metallicity Z cr is less than 10 -3 ' 5 Z Q , a value favored 
both theoretically (eg. Bromm et al. 2001, Schneider et al. 
2002) and observationally by the observed lack of metal- 
free stars in the Galaxy; (3) the spherical outflow model of 
Ostriker and McKee (1988). Thus it is fairly certain that 
PopIII formation continued well past the time at which the 
mean IGM metallicity reached Z CI , as large areas of the 
IGM remained pristine while others were enriched to many 
times this value. In fact, the history of PopIII formation is 
likely to have been determined almost exclusively by the 
distribution of high-redshift objects and the efficiency of 
metal ejection. 

To track the distribution of such high-redshift metals, 
we made use of the analytical formalism described in Scan- 
napieco & Barkana (2002) (for a similar approach see 
also Porciani et al. 1998). This formalism is a general- 
ization of the standard Bond et al. (1991) derivation of 
the halo mass function, and it interpolates smoothly be- 
tween all standard analytical limits, including the halo 
bias expression described by Mo & White (1996) and the 
Lacey & Cole (1993) progenitor distribution. By com- 
bining this general tool with an equally general charac- 
terization of the underlying IMF of primordial stars, we 
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were able to study the development of PopIII objects for a 
wide range of characteristic masses of these stars, ranging 
from 100 to 1000 Mq. This range was then further lim- 
ited by the simple observation that PopIII stars are not 
forming at z = 0. This is because choosing a characteristic 
mass in which a large majority of stars are outside of the 
140M Q < M* < 260M© mass range corresponding to the 
progenitors of SN 77 , results in very inefficient IGM en- 
richment, which allows PopIII star formation to continues 
indefinitely Thus lack of PopIII star formation in the lo- 
cal universe argues for a PopIII IMF with a characteristic 
mass of 200M Q ^ M c <; 5OOM . Even for the models 
in this range, however, the peak of PopIII star formation 
occurs relatively late (z ~ 10), and such PopIII stars con- 
tinue to contribute appreciable to the overall SFR well 
into the observable range of z <; 6. In all cases PopIII 

objects tend to be in the 10 6 ' 5 M Q — 10 7 M Q mass range, 
just large enough to cool within a Hubble time, but small 
enough that they are not clustered near areas of previous 
star formation. 

Employing a simple model of burst-mode star formation, 
we have developed a rough characterization of the directly 
observable PopIII objects. If PopIII stars formed accord- 
ing to a top-heavy IMF, the short but bright evolution of 
such stars would have boosted their luminosities by over 
an order of magnitude. Thus, PopIII objects would make 
up a large fraction of the detected high-rcdshift Lya emit- 
ters. In fact, for all top-heavy IMF models considered, 
we found that PopIII objects are potentially detectable at 
all redshifts beyond 5, although the precise fraction was 
sensitive to the assumed feedback efficiency. 

Such detections would be extremely difficult, however, 
if PopIII stars formed according to a normal IMF, for 
although metal-free stars are slightly brighter than their 
PopII/I counterparts (eg. Tumlinson, Shull, & Venkatesan 
2002) , this difference is negligible when compared to IMF 
effects. Thus in this case, PopIII objects would be de- 
tectable only if Eg 11 is small, although the absence of such 
galaxies can be used to set lower limits on the efficiency 
of feedback processes in the early universe. We note, how- 
ever, that the large equivalent widths seen in the observed 
populations of Lya emitters is well-explained by models 



containing a high fraction of PopIII stars (Malhotra & 
Rhoads 2001). If the IMF or feedback efficiency allows for 
direct detection, our model suggests the following obser- 
vational strategy: (1) collecting a large sample of high- 
rcdshift emitters; (2) selecting those at higher redshifts 
and closest to the favored flux range of 1.5 x 10~ 17 ergs 
cm~ 2 s _1 ; (3) selecting candidates based on high equiv- 
alent widths/lack of continuum detection; (4) performing 
deep, long-exposure spectroscopy aimed at the detection 
of the 1640A recombination line of Hell. 

Finally, our models allow us to asses the impact of 
PopIII objects on the ICM. We find that regardless of the 
assumed IMF, no more than 10% of the gas in clusters 
could have cooled into PopIII objects. Thus, contrary to 
previous studies (eg. Lowenstein 2001), we conclude that 
PopIII stars can contribute no more than 20% of the iron 
observed in galaxy clusters. However, if they formed with 
characteristic masses of about 200 — 260M Q , their pecu- 
liar elemental yields help to reconcile the observed Fe and 
Si/Fe abundances. Yet, even in this case, they overpro- 
duce sulfur and can account only for the O/Fe ratio in the 
inner regions of poor clusters. Furthermore, their energy 
input is insufficient to heat the ICM to the observed levels 
of 1 keV per particle. 

Early searches for PopIII stars have a long history of 
painful disappointment, driving astronomers to adopt a 
wide range of alternative techniques to explore their im- 
prints on later generations. Yet the same massive star 
formation that hid these objects from past searches may 
be the key to present day detections. We may finally be 
poised to make our first direct measurements of the pri- 
mordial generation of cosmic stars. 
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Fig. 1. — Cosmological metal dispersal. Top: Comoving radii of winds emanating from PopIII (solid lines) and PopII/I (dashed) objects. In 
the left panel we take Eg 11 = 10~ 4,0 and z a = 10 and consider sources at three mass scales, taken to be (from top to bottom) 10 9 Mq, 10 8 Mq, 
1O 7 M0, corresponding to fluctuations of initial sizes 190kpc, 90kpc, and 74 kpc, respectively. In the right panel we take £' n = 10" 2 5 and 
z s = 20 and consider 10 8 Mq, 1O 7 M0, and 1O 6 M0 objects, the smallest of which corresponds to fluctuations with a radius of 19 kpc. Bottom: 
Mean gas metallicity within each bubble. The lines represent the same models as in the upper panels, with Z w i n d increasing as a function of 
source mass. 
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Fig. 2. — Evolution of PopIII objects. Left: Fraction of PopIII objects for the best-guess models as a function of formation redshift for 
various masses and outflow energies. Each panel is labeled according to the energy input per gas mass, £g ■ In each panel the solid line, 
dotted, short dashed, long dashed lines indicate objects of masses 4.0 X 10 6 M Q , 1.0 X 10 7 M Q , 4.0 X 10 7 M o , and 1.0 X 10 8 Af Q respectively. 
The sharp drop in //jr at the smallest masses is due to the M m i n cutoff. Right: Fraction of PopIII objects for the conservative models, curves 
are the same as in the left panels. 
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Fig. 3. — Star formation rate densities for PopIII and PopII/I objects. Top Left : Star formation rate density per cubic comoving Mpc 
for PopIII objects (dashed lines) and PopII/I objects (solid lines). In this panel we assume /j = /j = 0.1 and consider models in which, 
from top to bottom, log(£g /z ) = —4.0 —3.5, —3.0, and —2.5. Bottom Left : Star formation rate densities in models in which = 0.5, for 

slightly larger log(5| 77 ) values of -3.5 though -2.0. Top Right: IGM average metallicity from PopIII (dashed lines) and PopII/I (solid lines) 
objects. Here f w = 0.3 and log(£'g // ) ranges from —4.0 to— 2.5 from bottom to top, as increasing feedback corresponds to additional SN 77 , 
which generate more metals. Bottom Right: IGM average metallicity with f w = 0.1, lines are as above. 
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Fig. 4. — Ionizing photon rate per unit stellar mass formed as a function of time in various (burst) models of star formation. The dotted 
line (labeled II) corresponds to PopII/I stars with metallicity of 0.05 solar and a Salpeter IMF for a range of stellar masses from 1 — 120Mq. 
The solid lines correspond to PopIII stars. The extended curve (labeled IIIN) is for a Salpeter IMF for a range of stars with masses from 
1 — 5OOA/0, while the peaked curve (IIIH) is for a "top heavy" Salpeter IMF, with a range of stars with masses from 50 — 5OOA/0. Finally, 
the dashed curve is the luminosity per solar mass for a single metal-free star with a mass of 260Mq . 
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Fig. 5. — Fraction of PopIII objects as a function of Lyo luminosity and redshift. Isocontours of fractions > 10 — 2 , 10 — 1 , 10 _1 and 10 ' 5 
are shown. Burst-mode star formation with a = fl 11 = 0.1 is assumed for all objects, and we adopt model II for Q(H) in PopII/I 
galaxies as given in Figure 4. In the PopIII case however, we assume a lower cutoff mass of 50Mq as given by model IIIH , in this figure. 
Each panel is labeled by the assumed Sg 11 value. For reference, the dashed line gives the luminosity corresponding to an observed flux of 
1.5 X 10~ 17 ergs cm -2 s — x , and the various points correspond to observed galaxies. The filled diamond is from Dey et al. (1998), the filled 
triangle is from Hu et al. (1999), the filled star is from Ellis et al. (2001), the open diamond is from Dawson et al. (2002), the open square is 
from Hu et al. (2002), the asterisks are from Lehnert & Bremer (2002), the open triangles are from Rhoads et al. (2002), the filled circles are 
from Fujita et al. (2003), and the filled squares are from Kodaira et al. (2003). Observations by Weyman (1998) and Ajiki (2002) are omitted 
for clarity of presentation. The curves have been extended slightly past z = 4 for comparison with the Fujita et al. (2003) data-set. 
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Fig. 6. — Fraction of PopIII objects as a function of Lya luminosity and redshift. In this case a Salpeter IMF with cutoff masses of 
1 — 5OOM0 is taken for PopIII objects, greatly lowering their luminosities. The dashed lines and points are as in Figure 5. 
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Fig. 7. — Number counts of Lya emitters per unit redshift and square degree as a function of redshift in various models of outflow generation. 
Top: Number counts assuming a top-heavy IMF as in Fig. 5. The dashed lines give the total number of Pop III emitters in models with 
assumed logio^g^) values (from bottom to top) of —2.5, —3.0, —3.5, and —4.0. The solid line gives the number counts of PopII/I objects 
in the log^^g 11 ) = —3.0 model, but is largely independent of this value. A threshold of 1.5 X 10 _17 ergs cm _2 s _1 is assumed in all cases. 
Bottom: Curves as in the upper panels, but assuming a Salpcter IMF as in Fig. 6. The filled circle gives the value observed by Hu et al. 
(1998). 
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Fig. 8. — Distribution of Lyo PopIII observed equivalent widths at two characteristic redshifts of 4.5 (upper panel) and 5.7 (lower panel) 
for a top-heavy IMF (solid), and two models assuming a 1 - 500M Q Salpeter IMF with log(£^. u ) = -3.5 (dotted) and log(f| /7 ) = -4.0 
(dashed). The upper panel in this figure is directly comparable to a possible PopIII contribution to the distribution in Fig 1. of Malhotra & 
Rhoads (2002) 
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Fig. 9. — Observed Fe abundance and O, Si and S abundance ratios relative to iron as a function of cluster richness. The data are 
normalized to solar photospheric values from Anders & Grevesse (1989). Each data point represents an average over clusters with similar 
richness. Triangles indicate ASCA observations and are taken from Table 2 of Fukazawa et al. (1998); dots indicate the results of a recent 
analysis of ASCA archival data by Baumgartner et al. (2002); the solid regions represent the average values observed by XMM (Tamura et 
al. 2002). 




Fig. 10. — Observed Fe abundance and O, Si and S abundance ratios relative to iron averaged over three cluster richness bins: [2-4] keV, 
[4-8] keV; [8-12] keV. The data are normalized to solar photospheric values from Anders & Grevesse (1989). 
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Fig. 11. — Top: Fraction of cluster gas mass that has collapsed into PopIII objects as a function of energy input; solid lines correspond to 
the best-estimate model, and the dashed lines correspond to the conservative model. In each model the lines indicate (from top to bottom) 
clusters with collapse redshifts of 2, 1.5, 1, and 0.5 respectively. Even in the lowest-energy, best-estimate model at most 10% of the cluster gas 
is processed by PopIII objects. Bottom: An optimistic estimate of the level of ICM heating due to PopIII outflows. Here we take f w = 0.1 
for all values of Eg 11 but assume that all of the supernova energy S^ II f^ 1 , is injected into the ICM. Even in this case the energy input in all 
models falls far short of the observed 1 kev/particle level of preheating. 
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Fig. 12.— Top Panels: Predicted ratio TZ = M^'/M* as a function of the energy per unit gas mass El assuming clusters formation 
redshifts 0.5 (lower curve) and 2 (higher curve). Horizontal lines indicate the predicted Eg 11 for models SN 77 -B (dashed), SN 77 -C (solid) 
and SN 77 -D (dot-dashed) with values for cr™ ln and <r™ ax as given in Table 2 assuming two different values for the PopIII star formation 
efficiencies: /» = 0.1 (left panel) and = 0.5 (right panel). The intersection with the predicted 1Z = Ml 11 /M* indicate the range of 

values for this parameter appropriate for each model. Bottom Panels: Predicted iron abundance in different SN 77 models (with the same 
line coding as in top panels) as a function of TZ. The solid regions indicate the range of Fe abundances that might be achieved in different 
models assuming the appropriate range of values for 1Z (see text). Left and right panels reflects the two different star formation efficiencies. 



24 




Fig. 13. — Predicted ICM (Fe/H), (Si/Fe), (O/Fe) and (S/Fe) abundances (in solar values) as a function of the relative number of SNIa 
and SNII. Dashed regions indicate the range of observed values in three cluster richness bins as in Fig. 10 and the different lines correspond 
to different SNII progenitor models (see Table 1). 
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Fig. 14. — Predicted ICM (Fe/H), (Si/Fe), (O/Fe) and (S/Fe) abundances (in solar values) as a function of the relative rate of Type la 
over Type II SNe assuming model SN 77 -B for PopIII stars and model SNII-B for Type II SNe. The long and short-dashed lines correspond 
to the limiting values for 1Z allowed by the model (see Fig. 12). For comparison, the solid line shows the predicted abundances if only Type 
la and Type II SNe contribute to the ICM enrichment. Dashed regions indicate the range of observed values in three cluster richness bins as 
in Fig. 10. 



